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Locality and stability of the cascades of two-dimensional turbulence 



oo 

o 

O 

(N 

(D 



Eleftherios Gkiouleka£| 

Department of Mathematics, University of Texas-Pan American , Edinburg, TX, United States 

We investigate and clarify the notion of locality as it pertains to the cascades of two-dimensional turbulence. 
The mathematical framework underlying our analysis is the infinite system of balance equations that govern the 
generalized unfused structure functions, first introduced by L'vov and Procaccia. As a point of departure we use 
a revised version of the system of hypotheses that was proposed by Frisch for three-dimensional turbulence. We 
show that both the enstrophy cascade and the inverse energy cascade are local in the sense of non-perturbative 
statistical locality. We also investigate the stability conditions for both cascades. We have shown that statisti- 
cal stability with respect to forcing applies unconditionally for the inverse energy cascade. For the enstrophy 
cascade, statistical stability requires large-scale dissipation and a vanishing downscale energy dissipation. A 
careful discussion of the subtle notion of locality is given at the end of the paper. 
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I. INTRODUCTION 

The physical notion of locality goes back to the 
Kolmogorov-Batchelor idea lUl 0. Ht] of an eddy cascade in 
three-dimensional turbulence where most of the energy is 
passed on from large eddies to smaller eddies by cascading 
through the intermediate scales. The dimensional analysis ar- 
gument behind the theory of two-dimensional turbulence pro- 
posed by Kraichnan |4] , Leith |53 and Batchelor ^ (KLB) is 
based in part on the conjecture that a similar physical princi- 
ple governs the upscale transfer of energy and the downscale 
transfer of enstrophy. In spite of the importance of the concept 
of locality to the foundations of the theory of hydrodynamic 
turbulence, there is no consensus on how to handle the con- 
cept rigorously. The need for a more rigorous understanding 
of locality becomes more pressing in light of some paradoxi- 
cal aspects of the theory of two-dimensional turbulence which 
will be briefly reviewed below. Becau se q uasi-geostrophic 
models of geophysical flows |0, HI HI [13. lUll relevant both to 
meteorology and oceonography, and two-dimensional mod- 
els of magnetically confined plasma turbulence lfl2i [Tsl [Till 
have a similar mathematical structure with two-dimensional 
turbulence, we cannot simply disregard the paradoxes of two- 
dimensional turbulence as irrelevant on the grounds that it is 
a fictitious fluid. 

For example, recent numerical simulations |[l5[ [l6l [TtI [isll 
have validated the KLB prediction for the energy spec- 
trum of the downscale enstrophy cascade. It remains unclear, 
however, whether the enstrophy cascade is a local cascade or 
nonlocal cascade. One side of the argument is that it cannot 
be a local cascade because the slope of the energy spectrum is 
too steep. On the other hand, if it is not a local cascade, then 
one has to explain why the prediction of dimensional analysis 
agrees with numerical simulations. Furthermore, it is worth 
remembering that prior to the groundbreaking paper by Lind- 
borg and Alvelius ITsll . every attempt to simulate an enstro- 
phy cascade failed. It is now understood that the presence of 
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a dissipation sink at large scales is necessMy for a successful 
simulation of the enstrophy cascade JT9I [20I [21I1 . Nonethe- 
less, we do not have a good grasp on why the presence of such 
a dissipation sink is sufficient. A recent theory by Falkovich 
and Lebedev ll22[ I23I1 predicts the scaling of the logarithmic 
corrections to the energy spectrum as well as the higher order 
structure functions of the vorticity for the enstrophy cascade. 
However, locality, and the existence of the enstrophy cascade 
itself are assumptions that are being entered into the theory. 
The relevant question is to understand theoretically the condi- 
tions needed for the existence of the enstrophy cascade. 

Ironically, the inverse energy cascade presents with an even 
more confusing situation. From a theoretical standpoint one 
would not expect the inverse energy cascade to be anything 
but local. From the standpoint of numerical simulations, there 
are m any pos itive re port s of the predicted k^^/^ energy spec- 
trum ||23[2l[2i[2l|2i[2l. The most convincing simulation 
of the inverse energy cascade has been reported in the paper 
by Boffetta et a/. ll29l1 , where in addition to the k^^/^ predic- 
tion, the 3/2 law has also been confirmed. On the other hand, 
the locality of the inverse energy cascade has been challenged 
on the grounds of numerical simulations giving conflicting re- 
sults OOl [31I [33, [33I1 . The current understanding is that un- 
der certain conditions there are coherent structures that spon- 
taneously form while the inverse energy cascade converges 
to stationarity. Apparently, the inverse energy cascade, as a 
physical process, continues to take place but it is hidden by 
the coherent structures which give the dominant contribution 
to the energy spectrum. Removing the coherent structures ar- 
tificially by postprocessing simulation data recovers the k^^^^ 
energy spectrum Isol [32iT34ll . This aspect of the inverse en- 
ergy cascade is not well understood. Furthermore, this phe- 
nomenon of the spontaneous generation of coherent structures 
is of considerable interest to oceonographers. 

In both cases reviewed above the issue at hand is the break- 
down of locality. The theoretical challenge is to understand 
how and why it happens. It should be noted that recent theo- 
retical work i2l[2l[3l[M[3l[3l[3l,[43] that expands on 
the KLB theory takes locality as well as the existence of the 
enstrophy cascade and the inverse energy cascade as assump- 
tions. As a result, although various aspects of these cascades 
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have been explained, the more fundamental question of the 
conditions needed for the existence of the cascades remains 
elusive. 

In the present paper we analyze the locality of the cascades 
of two-dimensional turbulence by adapting and gen eralizing 
the non-perturbative theory of L'vov et al. (EillH S, H 
Hslli^l- The mathematical framework is an infinite system of 
equations that govern the generalized unfused structure func- 
tions, the so-called balance equations. We also employ a 
scaling assumption, the fusion rules, which we conjecture to 
be valid in the enstrophy cascade and the inverse energy cas- 
cade. The fusion rules govern the scaling of the generalized 
structure functions when a subgroup of coordinates of veloc- 
ity differences approach each other. In previous work ll39il40ll . 
we used the balance equations to predict a linear superposi- 
tion principle between the downscale enstrophy cascade and 
the hidden downscale energy cascade which exists for finite 
Reynolds number. In that argument we did not use the fusion 
rules but we did assume the existence of the cascades. In the 
present paper we will consider more carefully the implications 
of the fusion rules on the existence question. 

The physical intuition behind our argument is as follows. 
Let Fn be the generalized structure function and let C„ be its 
scaling exponent. These structure functions satisfy a system 
of equations of the form 

Oji-Frj+i + /„ = DnFn + Qn- (1) 

Here, 0„F„+i is the nonlinear term that includes the effects 
of pressure and advection, /„ is a term associated with the 
sweeping interactions, Qn is the forcing term, and D„ is the 
dissipation operator. From the fusion rules it can be shown 
that the integrals in 0„i^„+i are local under the following 
conditions: for the downscale cascade UV locality requires 
C2 > and IR locality requires Cn+i < C2 + Cn-i; for the 
upscale cascade UV locality requires C„ — Cn-2 > and IR 
locality requires Cn+i > C2 + Cn-i- These conditions can 
be shown to be satisfied by the Holder inequalities. It follows 
that the interactions represented by 0„F„+i are local and also 
self-similar with scaling exponent Cn+i — 1- 

The implication of this argument is that the nonlinear inter- 
actions accounted for by the term 0„F„+i are local both for 
the enstrophy cascade and for the inverse energy cascade. This 
notion of locality is called statistical non-perturbative locality 
ll43ll . However, non-locality, in a different stronger sense, can 
arise from the forcing term Q„ - Although we may demand 
that the forcing spectrum be confined to a narrow interval of 
length scales, it does not follow that the forcing term Qn will 
force the balance equations only at those length scales. For the 
case of gaussian forcing, we show that the scaling exponent of 
Qn is Qn ~ Cn-2 +92 with ^2 = 2 for the downscale enstrophy 
cascade and q2 < for the inverse energy cascade. It follows 
that to have true locality we need qn — (Cn+i ^1) > in 
the downscale enstrophy cascade and g„ — (Cn+i — 1) < in 
the upscale energy cascade. These conditions are needed for 
the statistical stability of the cascades with respect to forcing 
perturbations. 

It should be noted that nonlocality via the forcing term Qn 
is only one of a number of possible scenarios for losing local- 



ity. The sweeping term /„ and the dissipation term DnFn can 
also destroy locality under certain conditions. A preliminary 
discussion of the sweeping term /„ was given in a previous pa- 
per ll47ll . and the dissipation term will be discussed in a future 
publication. Finally, it is also possible to lose locality through 
violation of the fusion rules. In that case, the term OnFn+i 
itself would not be local. In the present paper we will show 
that the UV locality of the term OnFn+i is very robust, even 
under violation of the fusion rules. However the same cannot 
be said for the IR locality. Our viewpoint then is to consider 
first the problems that can arise in the favorable case where 
the fusion rules are valid, before examining the validity of the 
fusion rules themselves in more depth. 

The argument of the present paper supports the conjecture 
of strong universality ll48ll for the direct energy cascade of 
three-dimensional turbulence and the inverse energy cascade 
of two-dimensional turbulence. However, it definitely rules 
out strong universality for the downscale enstrophy cascade. 
Because the argument relies on the hypothesis that the fusion 
rules hold for the downscale enstrophy cascade and the in- 
verse energy cascade, it is not completely rigorous. On the 
other hand, the hypothesis can be investigated by numeri- 
cal simulation. The p = 2 fusion rule, which is the es- 
sential one with respect to the locality argument, has been 
proven ll49ll50ll5lll52ll for the direct energy cascade of three- 
dimensional turbulence, and there is further support by exper- 
iments ll5l[54i|55L[56i[57il58ll. For the related problem of the 
passive scalar 1521, the fusion rules have been proved for all 
p (dda and have also been confirmed experimentally ll6ll l62ll . 
The problem of two-dimensional turbulence is similar enough 
to both problems to make the hypothesis plausible. 

From a more philosophical point of view, one can say that 
the scaling relations implied by the fusion rules are in fact 
a generalized definition of the physical concept of a "cas- 
cade". As has been pointed out previously ll43ll . from a phys- 
ical standpoint, the fusion rules mean that the large scales are 
correlated with the small scales in a very particular way where 
the self-similarity characteristics of the flow at the small scales 
"forget" the ongoing physical processes at the large scales 
(and vice versa for the inverse cascade) which leads to univer- 
sal scaling. The present argument then establishes the consis- 
tency between locality and the scale correlations needed for 
universality. The conditions needed for this consistency are 
necessary conditions for the existence of the cascades them- 
selves. 

The paper is organized as follows. Section 2 reviews the 
generalized balance equations of the generalized unfused cor- 
relation tensors, the emphasis being on distinguishing the 
sweeping interactions from the local interactions. Section 3 
introduces and motivates our revisions of the Frisch frame- 
work of hypotheses as the first step towards a theory of 
two-dimensional turbulence. The main idea is replacing the 
anomalous sink hypothesis with a universality hypothesis, 
which implies the fusion rules hypothesis. In section 4, we 
extend and generalize the locality proof of L'vov and Procac- 
cia ll43ll to the cascades of two-dimensional turbulence. In 
section 5, we then turn to the question of cascade stability, 
with respect to random gaussian forcing. We find that the 
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inverse energy cascade is stable, but that the enstrophy cas- 
cade is only borderline stable, with stability improving as the 
downscale energy flux is taken to zero. Section 6 discusses 
various subtleties that arise from our investigation regarding 
the concept of locality. Some technical matters are relegated 
to the appendices. 



II. THE GENERALIZED BALANCE EQUATIONS 

We now begin by reviewing the theory of the generalized 
balance equations. These equations were first derived by 
L'vov and Procaccia ||43ll and they are the foundation of pre- 
vious work ll39ll40ll47ll as well as this paper The two features 
of the balance equations that we would like to stress in this pa- 
per are the separation of the interaction term into local interac- 
tions and sweeping interactions, and the fact that the forcing 
term can be written in closed form for the case of Gaussian 
forcing. We also derive the balance equations that govern the 
generalized structure functions of the vorticity. 



A. Preliminaries 

The governing equations of two-dimensional turbulence 
are: 

+ UpdpUa = -daP + DUa + fa, (2) 
daUa = 0, (3) 

where is the forcing term, and D is the dissipation operator 
given by 



(4) 



Here the integers k and in describe the order of the dissipation 
mechanisms, and the numerical coefficients i/^ and /3 are the 
corresponding viscosities. D is the overall dissipation opera- 
tor The case k = 1 corresponds to standard molecular vis- 
cosity. The term fa represents stochastic forcing that injects 
energy into the system at a range of length scales in the neigh- 
borhood of the integral length scale £q. The term /3W^'^™Ua 
describes a dissipation mechanism that operates on large-scale 
motions. The operator V^^'" represents applying the inverse 
Laplacian repeatedly m times. In Fourier space this op- 
erator is diagonalized, and its definition may therefore be ex- 
tended to fractional values for m. The same holds for the 
parameter k. 

To eliminate pressure we multiply both sides of the Navier- 
Stokes equation with the operator Tap = Sap — dadpV^^ 
and we employ Tapup — up and Tapdp = to obtain 



dua 
~dt 



+ TapdjiupU^) = DUa + Tapfp- (5) 



Pa/3(x) as 

Tapvpix) = J dyPapiy^ - y)vp{y) 
= j dyPap{y)vp{x- y). 
For two-dimensional turbulence Pap (x) is given by 



Papi^) = SapS{x) - — 

Zn 



Oap 



Xa^P 



(6) 
(7) 

(8) 



The scalar vorticity ( is given by ^ = £apdaUp with Sap 
the Levi-Civita tensor in two dimensions. From the incom- 
pressibility condition daUa =0 it follows that there is a func- 
tion ip, called the streamfunction, such that Ua ~ Eapdptp. 
Using the identity Sap^p-y ~ Say one then shows that ( ~ 
Sap^p-ydadyijj — V'^ip from which we get ip — V^^C 

Ua = EapdpV^'^C- 

The vorticity equation is obtained by differentiating C, with 
respect to time and employing the Navier-Stokes equations: 



where C) is the Jacobian defined as 

J{A,B)=eap{dpA){daB), 



(9) 



(10) 



and g = Sapdafp is the forcing term. The nonlinear term 
J = J{ip, C) has been obtained by employing the following 
argument 

J = EapdaTp^ydsiUyUs) = Eapdaiu^d-yUp] (11) 

= UydyC + ieapdaUy){dyUp) (12) 
^UydyC = Jii'X)- (13) 

The term {eapdaUy){dyUi3) represents vortex stretching, but 
in two dimensions it can be shown that 

{eapdaUj){djUi3) = 0, (14) 

by direct substitution of the vector components. 

B. The balance equations 

To write equations concisely, we introduce the following 
notation to represent aggregates of position vectors 



X=(x,x'), 
{X}„ = {Xi, X2, . . . , X„}, 

{X}n = {^li • • • I Xfe_i, Xfc+i, . 



, x„}. 



(15) 
(16) 
(17) 



The operator Tap can be expressed in terms of a kernel 



We use the notation {X}„ + Ax as a shorthand to represent 
shifting all the constituent vectors of {X}„ by the same dis- 
placement Ax. Similarly, A{X}„ represents taking the scalar 
product of A with every vector in {X},i. Finally, the notation 
||{X}„|| ^ i? means that all point to point distances in the ge- 
ometry of velocity differences {X}„ have the same order of 
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magnitude R. And, the notation ||{X}„|| <C ||{Y}„|| means 
that all the point to point distances in {Y}„ are much larger 
than all the point to point distances in {X}„. 

Let Wa(x, x', t) be the Eulerian velocity differences 



The generalized balance equations can be derived by dif- 
ferentiating Fn with respect to t and substituting the Navier- 
Stokes equations (see appendixlAlfor details). This yields the 
equations 



Wa(x, x', t) = Uq(x, t) - Mq(x', t). 



(18) 



The eulerian one-time fully unfused correlation tensors are 
formed by multiplying n velocity differences Wc((x, x',t) 
evaluated at 2n distinct points 



OK 

dt 



OnFn+l + In = T)nFn + Qn- (20) 



Here !>„ is the differential operator representing dissipation, 
given by 



F„({X}„,t) = 



Lfe=l 



(19) 



When all velocity differences share one point in common, that 
is x'fe = Xo, we say that the correlation Fn is partially fused. 



J 



2)n - E['^(^x: + ^IV + /3(V.f" + V^,^")], (21) 

k=l 

and On is the linear integrodifferential operator such that 



(0„i^„+i)({x,x'}„,i) = J 0„({X}„,{Y}„+i)F„+i({Y}„+i,<) d{Y}n^ 

11 

= ^-DA;„({x,x'}„,t) = Dn{{x,X.'}n,t), 



(22) 
(23) 



fc=i 



where Dkn is given by 



1 " r 

i?,"r-""({x,x'}„,o = / dyp^MDzr''''''-''---'''^i{^y}n,y,t), 

1=1 



(24) 



with Dkni = Dknii + DknV2 + DknU + Dknii. and 



D 



knll 



({x,x'}„,y,i) = aa„+i,Xfc^;i+i' 



0!l---afc-l/3afc+i---Q 



D 



knl2 



({x, X }Ti,y,i) — 9Q,„_|_j_Xfc-F'„+i 



({Xm},„=i,Xfc -y,x'fc -y,{'KniYyn=k+i.'^k -y,X(), 

(25) 

({Xm}Jril>Xfc -y,x'fe -y,{l^m]\n=k+l^^k -y,x'i), 

(26) 



D 



«l---afc-l/9afc + i---a„ 



({x, x'}„, y, = a„„,,,x,,^^:;:"-^^"^+--""+^ ({X™},tii, X, - y, x', - y, {X™ 

(27) 



Ql---Qfc_l/3Q:fc + l---a 

A:ni4 



+--"({x,x'}„,y,t) = a„„,,,x,,F:- 



.Ql---afc_l^Qfc + l---Q„ + l 
+ l,X'fc J- „-f 



({X™},„Ji, Xfc - y, x'fc - y, {X™};\^fc+i, x'fc - y, x',) 

(28) 



The term /„ represents the sweeping interactions, and it is given by 

ja,a....a„ (1^^ ^^^^ ^ ^(a^,x. + a^,x' J (u,{{^k,^kU t) 

I 



fe=l 



i=l 



(29) 



where ^^({xfc, x'fc}„, t) is the generalized mean velocity; The term Qn represents the forcing term and it reads 



1 

U„({z,z'}„,t) = — EK(zfc,<) + u„(z'fc,t)), (30) 
" fc=i 



0„({X}„,t) = VQfe„({X}^,Xfe,t), (31) 



fe=i 
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where Qkn reads 



({X}„_i,Y,t) 



,fc=l 



(32) 

iP!3(Y,t)), (33) 



with(/j„(X,i) = fa(x,t) -/a(x',i). 



Let T„ be an abbreviation for the differential operator that 
transforms F„ to Vn such that Vn — TnFn- 

The balance equations for Vn and be derived easily by ap- 
plying the operator Tn on the balance equations for Fn- The 
result is 



C. Balance equations for the vorticity 

A similar set of equations can be derived for the general- 
ized structure functions of the vorticity. Let q(x, x', t) be the 
vorticity difference defined as 

g(x,x',t) = C(x,<)-C(x',t) (34) 
= ea/3(9Q,x + 9Q,x')w'/3(x,x',t), (35) 

and let y„({X}„,t) be the generalized structure function of 
the vorticity defined as 



K({X}„,0 = 



n9(Xfe,t) 



Lfc=l 



(36) 



It is easy to see that the vorticity generalized structure func- 
tions are related to the velocity generalized structure functions 
by 



K({X}„,i) (37) 

n 



k=l 



dVn 

dt 



' + 3n — '^nVn + Qr, 



(39) 



Here Q„ is the forcing term and 3„ is the sweeping term. The 
forcing term reads. 



Q„({X}„,t) = ^Q,„({X}^,X,,t), 



fe=i 



Qfe„({X}^,Y,i) 



.k=l 



,9(Y,t) 



(40) 
(41) 



To calculate the sweeping term we use (fT4l) to cancel the 
vortex tilting contributions. With a little bit of algebra we find 
(38) that 



J„({X}„,t) = J|[£a,ft(9a,,x, +ac^,x-)]X!(^7.Xfc + ^-y.x' J ( ({Xfe , x'fe }„ , t) 



i=i 



k=l 



J|u;ft(X(,x';,t) 



.1=1 



k=l 



Y[q{xi,x.'i,t) 



.1=1 



(42) 
(43) 



The trick is to apply the operators ea^Pj (t^a,,x, + da^.^'j) 
one by one onto the ensemble average in eq. ( l42b . wherein 
ri — 1 of the factors are constant for j ^ I with respect to 
Xj, x'; , and use the identity £q.^9q.9^[m^m^] = u^d^C on the 
wp^ and lij factors that are both x^ , x'j dependent. Each ap- 
plication of these operators effectively converts each wp^ fac- 
tor into a corresponding ^(x;, x'j, t) factor The exact math- 
ematical form of the term TnOnTn+i^n+i is not required. It 
is only sufficient to note that once it is shown that the expres- 
sion OnTf^^iKi+i = 0„F„+i is local, then it easily follows 
that the term TnOn'^n+i^n+i is also local since T„ is a hnear 



differential operator. 



III. THE THEORETICAL FRAMEWORK 

Both the K41 theory for three-dimensional turbulence, and 
the KLB theory for two-dimensional turbulence are based on 
a dimensional analysis argument. However, Frisch ll63i l64ll 
has suggested that Kolmogorov's second paper 101 leads to 
the following more rigorous reformulation of the dimensional 
analysis argument, based on the following three hypotheses: 
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(HI): At small scales and away from any boundaries, the 
velocity field is incrementally homogeneous and incremen- 
tally isotropic; (H2): Under the same conditions, the veloc- 
ity field is self-similar at small scales, thereby possessing a 
unique scaling exponent h ; (H3): the turbulent flow has a 
non-vanishing mean dissipation rate in the limit of infinite 
Reynolds number (i.e., an anomalous energy sink). Then, one 
uses (HI) and (H3) to derive the 4/5 law which implies that 
/i = 1/3, and from (H2) the scaling for all structure functions 
and the energy spectrum is deduced. 

In a recent paper, Frisch ll65ll questioned the self- 
consistency of the assumption of local and incremental ho- 
mogeneity. The argument essentially is that it is not obvious 
whether the nonlinearity of the Navier-Stokes equations will 
preserve incremental homogeneity unless the initial condition 
is globally homogeneous. In a previous paper ll47ll I have ar- 
gued that incremental homogeneity will be preserved in the 
upscale and downscale inertial ranges only if the sweeping 
interactions, represented by the /„ term of the balance equa- 
tions, can be neglected in the inertial range. As I have empha- 
sized in that paper, this condition on the /„ term is necessary 
for the very existence of an inertial range! Here we will sim- 
ply take it for granted in order to focus our attention on the 
other needed conditions. 

Within the Frisch framework, many theoretical approaches 
to three-dimensional turbulence that try to predict the inter- 
mittency corrections to the scaling exponents of the structure 
functions, can be interpreted as extensions of the Frisch theory 
where the self-similarity assumption (H2) is weakened while 
the other two assumptions (HI) and (H3) are tolerated. It is 
an easy exercise to reformulate the dimensional analysis ar- 
gument of the KLB theory in a similar manner. However, a 
theory along these lines would already take for granted the lo- 
cality and universality of the two cascades. Contrary to the 
situation in three-dimensional turbulence, what we must un- 
derstand are the conditions needed to satisfy universality and 
locality. In previous work we have proposed that the 

questions of locality and universality can be probed more rig- 
orously by adapting the theoretical work of L'vov and Pro- 
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caccia et al. ||41l |42|, |41 |4J, |4S IM |52|] to two-dimensional 
turbulence. We will now expand further on this idea on the 
remainder of the present paper. 



A. Revisions to the Friscli framework 

We propose that the Frich framework of hypotheses should 
be revised as follows: 

First, we adopt Frisch's (HI) to our formulation. We have 
shown previously |4'!^i| that a stronger homogeneity hypothesis 
is needed to eliminate the sweeping interactions before deriv- 
ing the 4/5 law. Though we may circumvent this problem by 
postulating that stronger assumption of homogeneity for our 
framework, we believe that it is desirable to be able to estab- 
lish the stronger hypothesis from first principles (see section 
5ofRef.||47ll). 

Second, to allow for intermittency corrections, it is neces- 
sary to relax the self-similarity hypothesis (H2). One pos- 
sibility is the space-time self-similarity assumption, used in 
the early papers of the quasi-Lagrangian diagrammatic the- 
ory ll66ll67ll . It was shown later that this assumption is false, 
because it axiomatically implies Kolmogorov scaling and for- 
bids intermittency corrections ll44ll68ll . thus leading to a self- 
inconsistent theory. The successful proposal is space one-time 
self-similarity, defined below, and we shall adopt it in this pa- 
per 

Third, following L'vov and Procaccia ll4lll42ll43ll . we adopt 
an hypothesis of universality. Its purpose is to replace the ad 
hoc assumption of anomalous sinks. The universality hypoth- 
esis itself claims that statistical symmetries are recovered at 
length scales away from the forcing range even when the en- 
semble is constrained by a symmetry-breaking condition at 
scales closer to the forcing scale. 

Taking the ideas above into consideration, we postulate the 
following hypotheses for both the enstrophy and energy iner- 
tial ranges: 



Hypothesis 1: The velocity field is incrementally stationary, incrementally homogeneous, and incrementally isotropic, defined 



afi,({X}„,t) 
dt 



= o,vt e 



fc=i 



F„({X}„,t) = ^^„(ro +yi({X}„ - ro),<), G 50(2). 

as long as the evaluations {X}„, {X}„ + Ar, Tq + yi({X}„ — Tq), lie within an inertial range. 

Hypothesis 2: The velocity field is self-similar in the sense that for every evaluation {X}„ within an inertial range 

Be > : F„(A{X}„, t) = A^"F„({X}„, t), VA e (1 - e, 1 + e). 

For the hypothesis of universality, we define the conditional correlations 



$„({X}„,{Y},„,{w,.}^Li,i) 



Lfc=l 



(44) 
(45) 
(46) 

(47) 
(48) 



7 



and use them to formulate the additional hypothesis that in the inertial range, the conditional correlations $„ essentially honor 
the same symmetries as the unconditional correlations Fn, in the asymptotic limit where ||{Y},„|| are situated between ||{X}„|| 
and the forcing scale Iq. 

Hypothesis 3: Let {X}„ and {Y}„i represent the geometries of velocity differences and let W = W({Y}m, {'^kYk=i)- 
Then, if in the direct cascade they satisfy ||{X}„|| <C ||{Y}m|| <C ^o. or alternatively if in the inverse cascade they satisfy 
||{X}„|| ^ ||{Y}„i|[ ^ £o. then the conditional correlations preserve incremental stationarity, incremental homogeneity, 
and incremental isotropy, with respect to {X}„, defined as 

n 

E(^".,x. +5a.,x'J<f>„({X}„,{Y}„„{w,}r.i,i) = (50) 

k=l 

4>„({X}„,{Y}„„{wfeK!Li,t) = $„(ro +yL({X}„ - ro), {Y}™, {w,.}^Li, i), V^L e 50(2), (51) 
and also self- similarity, with the same scaling exponents (^„, defined as 

3£ > : <1>„(A{X}„, {Y}™, {wfel^Li, i) = Af"$„({X}„, {Y},„, {wfc},"Li, i), VA e (1 - e, 1 + e). (52) 



Hypothesis 1 is essentially the first hypothesis in the 
Frisch formulation. Hypothesis 2 is the space one-time self- 
similarity principle introduced by L'vov and Procaccia ||43I1 
in the context of three-dimensional turbulence. The scaling 
exponents represent the scaling structure of each inertial 
range. If < C2 < 2, then the energy spectrum follows a 
power law given by E{k) ~ k~^~''^ ll64ll . If there is a loga- 
rithmic correction, then the result also holds for <^2 = 2. Hy- 
pothesis 3 states that the statistics of the velocity field at a 
certain scale still maintain the symmetries stated in hypothe- 
ses 1 and 2 even when a symmetry-violating constraint is im- 
posed via a conditional average at scales closer to the forc- 
ing scale. The constituent statements of hypothesis 3 shall be 
referred to as universal incremental homogeneity, universal 
incremental isotropy, and universal self-similarity. Note that 
it is essentially a more careful reformulation of the assump- 
tion of "weak universality" that was proposed previously by 
L'vov and Procaccia ll4lll43ll . The underlying idea is that the 
condition w{xk,x'k, t) = in the definition of the condi- 
tional correlations (1>„ partitions the ensemble of all possible 
forcing histories consistent with the overall forcing spectrum 
and the stationarity assumption into subensembles defined by 
the parameters {'WkY^^i- Each choice of {Y},„ represents a 
distinct partition of the entire ensemble into subensembles. 
The assumption for the statistical behavior of the velocity 
field is that it remains invariant accross each subensemble of 
forcing histories for all subensemble partitions {Y}„i (with 
||{X}„|| <C ||{Y}„ii| ^ £0 if it is a downscale cascade or 
j|{X}„|| ':$> ||{Y}m|| ^ £0 if it is an upscale cascade), and 
thus dependent only on the overall forcing spectrum. 

B. The fusion rules hypothesis 

The immediate consequence of the universality hypothesis 
is the fusion rules, whose physical interpretation is that differ- 
ent length scales are correlated (a hint of the cascade process) 
and that the governing interactions, as we shall show in the 



next section, are local (a consequence of the fusion rules and 
the structure of the Navier-Stokes equations). In a forthcom- 
ing paper, we will show that the fusion rules also govern the 
location of the dissipation length scales and that, in doing so, 
they provide anomalous energy and enstrophy sinks! 

Consider a geometry of velocity differences {X}„ such that 
all point to point distances have order of magnitude 1, and 
define 

FiP)(r,i?) = F„(r{X,}Li,i?{Xfc}^=p+i). (53) 

The function F,f^ (r, R) reflects the case where p veloc- 
ity differences have separations with order of magnitude r, 
and n — p velocity differences have separations with order 
of magnitude R. The case of interest is when the evalua- 
tion ir{Xk}l^^,R{Xky^^ 

p+i) is within the inertial range 
3n Q (M^)^" and r <^ R. The fusion rules give the scaling 
properties of F^^^ in terms of the following general form: 

Fif)(Air, X2R) = xl"'Xi--^-''F^[P\r, R). (54) 

Since Fn is defined as the product of velocity differences we 
expect the limits Ai and A2 ^ to converge. This 
implies that > and C„ — ^„p > . A concise state- 
ment of the fusion rules hypothesis is that for the direct en- 
strophy cascade ^,ip = , and for the inverse energy cascade 
£,np ~ Cn — Cn-p for 1 < p < 71 — 1 . The cases p ~ 1 
and p ^ n — 1 require some additional considerations, and 
can be deduced, as it turns out, from the p ~ 2 fusion rule 
(see section lTV Ab . We will also consider the case of "regular" 
violations to the fusion rules where the scaling exponents 
satisfy < ^„p < Cn, so that the exponents on Ai and A2 are 
both positive. 

We will now briefly review the argument of L'vov and Pro- 
caccia |[43tl that that the fusion rules hypothesis is an imme- 
diate consequence of the universality hypothesis. Let us con- 
sider first the case of the direct enstrophy cascade. For the case 
2<p<n-2we will show that for ||{X}„|| < ||{Y}„j| the 
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fusion scaling is Let J'({X}„, {wfej^^j^) be the probability for the event 

w(xfc, x'fc, t) = Wfc. It follows that 

(55) 



F„(A{X}p,/x{Y}„_p) = J 



'n—p 
.k=l 



n—p 



T(^{Y}„_p, {^k}lZl)%{X{^h, /i{Y}„_p, {wfcllCf ) n dwk (56) 



n—p 
W Wk 

Lfc=i 



fc=i 



n^^{YU-p, {wfcl^Zf )$p({X}p, m{Y}„-p, {wfcl^Zf ) n (57) 



fc=i 



n—p J • 



(58) 



The factor Fn{{'K}p, ^{Y}„_p) is now independent of A and has to scale as /^''""^p. 

For the case of the inverse energy cascade, again for 2 < p < n — 2 and under the same limit ||{X}„|| ^ ||{Y}„|| the fusion 
scaling is 



i^„(A{X}p,Ai{Y}„_p) - A«"-^"-/--^^„({X}p,{Y}„.p). 
We show this with a similar argument as follows; 



K(A{X}p,^{Y}„_p) = J 



p 

.k=l 



nXiXjp, {^kYU)<i>„_p{f,{YU.p, A{X}p, {w4^Zf ) n dwk 



k=l 



P 
.k=l 



nH^}p, {w,.}L,)$„_p({Y}„_p, \{X}p, {wfc}[!Zf ) n dwk 
/"-F„(A{X}p,{Y}p). 



fe=i 



(59) 

(60) 

(61) 
(62) 



The factor F„ (A{X}p, { Y}„_p) is now independent of yU and 
has to scale as A^"~^"-p. 

We would like now to briefly discuss the motivation behind 
our conjecture that the enstrophy cascade and the inverse en- 
ergy cascade satisfy the fusion rules. First, it should be noted 
that for the locality proof given in section|IV]we only need the 
fusion rule for the cases p = 2 and p ^ n — 2, from which 
one then derives the scaling for the cases defined in Fig.[T]and 
Fig. 12] For the energy cascade of three-dimensional turbu- 
lence the p = 2 fusion rule has been demonstrated by Feyn- 
man diagram analysis ||49[ Isol Isil l52ll . The proof indicates 
that the fusion rule essentially follows from the assumption 
that the scaling exponent ^2 is universal and does not change 
in response to perturbations to the forcing statistics. This as- 
sumption rests on less solid ground for the enstrophy cascade, 
however we can expect it to be true at least in the experimental 
situations where the cascade actually exists. It is also worth 
noting that this assumption is weaker than our hypothesis of 
universality, which in some regard is a stronger assumption 
than what is really needed. 

There is another consideration that strongly motivates our 
conjecture: the p = 2 fusion rule controls the positioning of 
the dissipation length scale ||42[ l43i l54ll . In a forthcoming pa- 
per we will show that if this fusion rule is violated, then the 
dissipation length scale would not be correctly positioned to 



dissipate the injected energy or enstrophy. Consequently, it is 
not easy to reconcile the numerical observation of both cas- 
cades with a violation of the fusion rule p ~ 2. Furthermore, 
a situation where the p = 2 rule is satisfied and the other rules 
are violated is unlikely. Finally, in two-dimensional turbu- 
lence, due to the smaller dimensionality of the problem, we 
are afforded the opportunity to test of validity of the fusion 
rules directly with a numerical simulation. 



C. Symmetries and the balance equations 

The assumptions that we have put forth are not self-evident 
axioms but hypotheses. Thus, the goal of theory is not only 
to derive conclusions from these assumptions but to also work 
in the opposite direction and give reasons that justify the as- 
sumptions themselves. 

The argument that was given by Frisch ll63il64ll begins with 
the observation that the unforced Navier-Stokes equations are 
invariant with respect to space and time shifts and rotations: 

(f,x,u) ^ (t,x + Ax,u), VxeM'^, (63) 
(i,x,u) ^ (t,^x, Au),VA e S'O(d), (64) 
(i,x,u) -> + At,x,u). (65) 
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Furthermore, if we ignore the dissipation terms, then the 
Navier-Stokes equations are also invariant with respect to the 
following self-similar transformation 

(t,x,u) ^ {X^-''t,Xx,X''u)yX eR+,he R. (66) 

In hydrodynamic turbulence these symmetries are obviously 
broken by the forcing term, the boundary conditions, and 
the self-similarity symmetry by the dissipation terms. Frisch 
ll63[ l64ll hypothesized that these symmetries will be statisti- 
cally reinstated in the inertial range when the flow is governed 
by a strange attractor. The big question is; how do we prove 
this? We believe that the generalized balance equations, de- 
rived in the previous section, are the proper theoretical frame- 
work within which this question can be addressed. 

We begin by accepting the assumption of local stationarity 
for the reasons given by Frisch ll63l l64ll . Then, the balance 
equations read 

0„K+i + /„ = T>„F^ + Q„. (67) 

As was pointed out previously ll44[ l69ll . the advantage of us- 
ing generalized structure functions where every velocity dif- 
ference is associated with two distinct coordinates that are dif- 
ferent from any other velocity difference, is that in the limits 
1^ ^ and /3 — > the dissipation terms can be dropped. 
This is not possible for the standard structure functions where 
every velocity difference is associated with the same two co- 
ordinates. We show this by using the mean-value theorem to 
bound T>nF„ as follows 

|DnF„| < + C2mlZ?j \Fn\. (68) 

V ^min / 

Here, Ci and C2 are constants independent of ly and (3, and 

i?min = min{xfc, x'fe : k e N,l < k < n}, (69) 
i?niax = max{x/c, x'/c : /c S N, 1 < < n} . (70) 

It is easy to see that because all the differentiations can be 
performed without invoking the product rule, the viscosities 
1/ and [3 multiply on a factor that remains finite in the limits 
1/ — > and /3 — > 0. Thus, jl'n.F'nl in the inertial range. 

To reinstate the statistical symmetries we need a region of 
length scales where Qn and /„ can also be ignored. Then, one 
has the homogeneous equation OnFn+i = 0, which remains 
invariant both under local homogeneity and local isotropy. In 
fact, it is also known jiil l45ll that the homogeneous equation 
is invariant under the following similarity transformation: 

{X}„^A{X}„, ^ A""+^('')f„. (71) 

consequently, it is expected to have solutions in the general 
form 

J dfiih)F„M. (72) 

where Fn.h are the zero-modes of the operator 0„ which scale 
as 

K,,,(A{X}„,t) = A""+^('''F„,„({X}„,t). (73) 



Note that the same result can also be obtained from the multi- 
fractal hypothesis ll64ll . On a large inertial range, the leading 
contribution to Fn is asymptotically self-similar with the scal- 
ing exponent (n given by 

Cn = mm{nh + Z{h)). (74) 

h 

For the case of a multifractal stochastic velocity field with 
D{h) defined as the fractal dimension of the set of points that 
support a local Holder exponent h, the relationship between 
Z{h) and D{h) is Z{h) = d — D{h) where d is the dimension 
of the velocity field, and d = 2 for two-dimensional turbu- 
lence. 

It has been suggested, for the case of three-dimensional 
turbulence, that the scaling exponents d can be calculated 
from the solvability condition of the homogeneous equation 
0„F„+i = HlSllMl- Although, from a practical stand- 
point, perturbative methods have been more effective lf7ll[72ll . 
the solvability condition analysis reveals the underlying prin- 
ciple governing the origin of the scaling exponents C„. From 
a physical standpoint, the condition 0„i^„+i = includes 
(for n = 2) and extends (for 7i > 2) the requirement of a 
"constant" (in the asymptotic sense) energy flux in the inertial 
range. The extension makes the condition powerful enough to 
lock down all the scaling exponents (n, as was demonstrated 
by Belinicher et al. |7(||. As we have shown in a previous 
paper ll39ll . the problem with extending this argument to two- 
dimensional turbulence is that the scaling exponents C,n of the 
enstrophy cascade are not non-trivial solutions to the equation 
OnFn+i = 0. This should not surprise us, that we cannot 
obtain the scaling exponents of the enstrophy cascade from 
an "extended" constant energy flux condition! What must be 
done instead is to use the equation 

T„0„F„+i = 0, (75) 

obtained by the generalized balance equations for the vor- 
ticity structure functions derived previously in section III CI 
This equation represents an "extended" constant enstrophy 
flux condition, and it yields two solutions for the scaling ex- 
ponents, instead of just one: an energy cascade solution that 
transfers energy but not enstrophy (the non-trivial solution 
of OnFn+i = and it also satisfies T„0„i^„_|-i = triv- 
ially because it does not transfer enstrophy), and an enstrophy 
cascade solution that transfers enstrophy but not energy (the 
non-trivial solution of T„0„i^„_(.i = and it also satisfies 
0„F„_|-i = trivially because it transfers no energy). It also 
follows from the mathematical structure of the equation ( f75b 
that these two solutions can be superimposed linearly to obtain 
a composite solution that transfers both energy and enstrophy. 
The possibility and implications of such a c omp osite solution 
has been discussed in previous papers iH, Hfl, S and 
will not concern us further in this paper 

These observations show that a constructive point of view is 
to see our hypotheses 1,2, and 3 as an efficient definition of the 
concept of an "inertial range", in a generalized sense. Obvi- 
ously, the hypotheses are valid only on a multidimensional do- 
main of velocity differences geometries {X},i S tJ„. The ex- 
tent of this domain 3„ is the extent of the inertial range itself. 
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Xi = X2 



FIG. 1 : The p = 1 fusion rule geometry with a type B fusion. Here 
we take the limit r <^ R with r and R both in the inertial range. In 
a type B fusion, the small velicity difference shares an endpoint with 
one of the large ones, i.e. Xi — X2. 



R 

H 




FIG. 2: The p — n — 1 fusion rule geometry with a type B fusion. 
This is a composite rule where we take the limits I <^ R and r <C _R. 
The velocity difference associated with i shares an endpoint with the 
velocity difference associated with R. 



A one-dimensional interval of length scales where the struc- 
ture functions exhibit power law scaling, is a reduction of the 
domain tf„ in which information is lost. For the case of dual 
cascade, we have an upscale range and a downscale range, and 
a different set of scaling exponents and region Sn is asso- 
ciated with each range. To determine the extent of the region 
3n for the energy and enstrophy ranges we employ the theory 
of the generalized balance equations, combined with the fu- 
sion rules hypothesis. More rigorously, the domain 2n is the 
range of length scales where the terms Qn, In, and DnFn in 
the generalized balance equation are negligible relative to the 
terms contributing to 0„F„+i. The first step towards deter- 
mining the extent of the domain 3„ is to calculate, from our 
hypotheses, the scaling exponents of the terms of the balance 
equations. Then these terms can be compared against each 
other. We initiate this study in the next two sections of this 
paper. Note that it is sufficient to study in this manner only the 
balance equations for the velocity field. Since the operator T„ 
is a strictly differential operator, it is also local, therefore the 
scaling exponents of the terms of the vorticity balance equa- 
tions are equal to the the scaling exponents of the terms of the 



velocity balance equations minus 1 . So, the scaling exponents 
for pairwise ratios of the terns against each other are the same 
for both balance equations. 

IV. LOCALITY OF THE INTERACTION TERM 



We will now show that the p = 2 fusion rule and p = n — 2 
fusion rule combined with local homogeneity and incom- 
pressibiUty, implies that the nonlinear interactions in the in- 
ertial range are local. From the viewpoint of the generalized 
balance equations, the nonlinear interactions are accounted for 
by the integral in the term 0„F„+i, and the sweeping inter- 
actions by the term /„, which we assume,_for now, that it is 
negligible in the inertial range (see Ref. 47 for further dis- 
cussion). We say that the integral is local if it is convergent 
and furthermore if the dominant contribution to the integrals 
in OFn+i comes from the region in which the separation of 
the integral variable y from all other points has the same or- 
der of magnitude as {X}„ ^ R. Locality implies that the 
contributions Dkn to OF„+i are also self-similar with scaling 
exponent (5„ and satisfy 

-Dfc„(A{X}„,t) ^ A'5"i?fc„({X}„,t), (76) 

where 5„ is given by (5„ = Cn+i — I- We propose that the local- 
ity of the interaction integral in Dkn is the mathematical defi- 
nition that corresponds most closely to the kind of locality that 
is required to enable an eddy cascade with universal scaling. 
In the sense of our proposed definition, we will show that both 
the energy and enstrophy cascade of two-dimensional turbu- 
lence are local. 

The proof given in this section is based on a previous proof 
by L'vov and Procaccia given in section IV-C of Ref. |43l The 
same argument is also presented in the appendix of Ref. |6§. 
We have generalized their proof in two directions: first, we 
derive the explicit conditions needed for locality even for the 
case where the fusion rules do not hold; second, we extend the 
proof to the case of the inverse energy cascade. 

A. Preliminaries 

It can be seen from the equations (|25ll, (|26ll, (|27l), ( l28b 
that the general form of the terms that contribute to Dkn in- 
volves an integral of the form 



3 y" cJy f'«,/3(y)57,xfc 



l=l,l^k 



W0{xk - y,x'fc - y)w^{xk - y,s) ) , 



(77) 



r 



where s can be any point among xi, . . . , x„ or x'l, . . . , x'„. 0, x^ — y ^ x; or x';, x'^ — y ^ x/ or x'/, and p = ||y| 
The locality proof requires the scaling of F„ in the limits y 
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CO. Consequently, we need the fusion rules for the geometries 
shown in Fig. [T](case p ~ 1) and Fig. |2](case p ~ n — 1). 
Both can be derived from the fusion rules for the cases p = 2 
and p ^ n — 2. 

(a) For the case p = 1 where we also assume a type IB 
fusion (i,e. xi — X2, and see Fig. [TJ the governing fusion rule 
is 

Fn - {r/Rf^R'^^ (downscale), (78) 
F„ - (r/i?)'^"-^'-^i?«" (upscale). (79) 

To show this, we note that 

W(X2,X'2) = W(x2,x'i) +W(x'i,x'2) (80) 

= w(xi,x'i)+w(x'i,x'2). (81) 

For the last step, we used xi = X2. Let Y = (x'l, x'2). Then 

K({X}„) = 7^„(Xi, X2, {X};!^3) (82) 
= F„(Xi, Xi, {X}^^3) + ^^„(Xi, Y, {X};!^3). 

(83) 

The third term is the same fusion problem as the first term 
because Xi and Y share the point x'l, and from the universal 
isotropy hypothesis we can rotate the legs r and R in Fig. 
[T] with respect to each other so that the three points form an 
isosceles triangle. Then one problem can be obtained from 
the other problem by reflection around the triangle's axis of 
symmetry. Consequently, both problems scale according to 
the second term, which is a p = 2 fusion. In the proof below, 
we will use the generalized scaling 

F„ ^ (r/i?)«"-^i??", (84) 

which is applicable both upscale and downscale. 

(b) For p ~ n — \ with type B fusion, we have n — 2 ve- 
locity differences of order r, one velocity difference of order 
I with one endpoint attached to a velocity difference of or- 
der R, where £ <^ R and r ^ R. Note that this fusion 
can be composed as follows. Begin with all velocity differ- 
ences at order R. Then take the following limits: (ii) Shrink 
one velocity difference to order £ <C i? with one endpoint 
attached to another velocity difference (this is the previous 
case); {£2) Shrink all other ?i — 2 velocity differences down to 
order r <C R- Thus, we have, for the downscale case. 



The first limit (^i) gives the first factor (^/i?)^^, and the sec- 
ond limit (£2) the second factor (r/i?)^"-^ Similarly, for the 
upscale case, using the exact same hmits (^1) and (£2) , we 
find 

'""^ (^)'""%^- (87) 
~ £C"-C"-2/"-C2^C2+a-2-C„^ (88) 

In the proof below, we will use the generaUzed scaling 

B. UV locaUty 

UV locality requires convergence in the limits y 0, x^ — 
y ^ X; orx';, and x'^. — y ^ x; orx';. The only limit 
that requires serious consideration is the first where Pap (y) 
is singular For this case we distinguish the following two 
subcases. 

(a) Assume that x.^ ^ s. The derivative of the ensemble 
average in ( ITTI i is analytic in y ^ 0, so we Taylor expand it 
around y = 0. 



3- / dy Pa,fi{y)[Ap+'Bi3^y^ + ei3^syjys + ---]- (90) 



The first term vanishes by incompressibility. The second term 
vanishes because the integral is odd with respect to y, from the 
local isotropy hypothesis, whereas PakPiy) is even. The third 
integral is local. Use dy = p dp dQ,{A) with p = ||y||, A S 
50(2), and dVL{A) the measure of two-dimensional spherical 
integration. The third integral then reads 



h = j dp j dn{A)pP^,p{y)Qp^sy^ys (91) 
~ / dppp-'^p' f dpp^p\ (92) 

and it is unconditionally local 

(b) Assume x^ = s. Then the integral reads 



/ 



dy Pa^0iy)d. 



n 



wpi^k - y,x'fe - y)u;^(xA; - y,Xfc) 



(93) 
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and in the limit y ^ we have the velocity difference ge- 
ometry shown in Fig. [3] From the p = 1 fusion rule with 
type IB fusion, the ensemble average in the integral scales as 
Fn+i ~ The integral then scales as 



dp pp^V^V^"^'-" 



0+ 



dp p' 



,«« + i,2-2 



(94) 



Here, the spherical integral contributes the factor p, the pro- 
jection operator Pa^,/3 (y) contributes p~^, the derivative d^^^k 
contributes p^^ (because the x dependent factor depends only 
on the smallest in separation of the two velocity differences 
in Fig. [3] which makes that factor dependent only on p), and 
the fusion rule contributes p5"+i.2. The resulting integral is 
marginally local for £,n+i.2 = C2 ~ 2 (enstrophy cascade) and 
non-local for 12 = C2 = 2/3 (downscale energy cascade 
in 3D). However, note that the type IB fusion rule for the case 
p = 1, which we have used here, is written in more detail as: 

F„+i ~ (w/3(xfc - y,XA;)?i;^(xfe - y,Xfc)) (95) 
~ $2(xfc -y,Xfc,Xfe -y,Xfe)$„_i, (96) 

which allows the integral 3 to be rewritten as 



J dy PQ^/3(y)9^.Xfc$2(xfc - y,Xfc,Xfc - y,Xfc). 

(97) 

Here we have used the fact that is independent of both 

(xfc-y, Xfc) and (xfc-y, x'fc-y), thus independent of x^, and 
therefore it can be pulled out of the 9^,xfc operator It is easy 
to see that the leading term of the $2 factor vanishes when 
differentiated by 5^,Xfc by universal incremental homogeneity. 
Thus, we get a cancellation that kills the leading contribution 
and the integral then scales according to the next-order term: 



dp p' 



ri5" + 1.2 



(98) 



0+ 



This integral is local if > (i.e. for locality we 

need 3 ^ as p — > 0). The result holds unconditionally, 
even under a regular violation of the p = 2 fusion rule, e.g. 
Fn+i ~ as long as ^„+i,2 > and some 

factorization Fn+i ^ ^2^n-i is still possible (that would be 
true for higher-order terms, if the leading term should happen 
to vanish) . Under the fusion rules hypothesis this condition is 
C2 > for a downscale cascade and Cn+i — Cn-i > 0, Vn G 
N — {0, 1} for an upscale cascade. 

Consider finally the cases x^ — y x; or x';, and 
x'fe — y ^ X; orx';. We perform the integral spherically 
around the value of y where one of these coincidences take 
place. Let p be the distance between the two approaching 
points. Assume any regular fusion rule of the form Fn+i ~ 
Now, the function Pakpiy) is no longer 
singular so we gain a factor of p^. Otherwise, the computa- 
tion is the same as in the previous case, and the integral scales 
as 



dpp' 



+ + 1 _ 2+2 



(99) 



0+ 



which is local even under a regular violation of the p = 2 
fusion rule. 



C. IR locality 

Consider the limit p ~ ||y|| 00. The corresponding 
geometry of velocity differences is shown in Fig. |4l For the 
downscale cascade we use the fusion rule for the case p = 
n — 1, defined in Fig. |2l 



n+l 



Cn + 1,2 



R 



(100) 



Expanding around the point at infinity p ^ 00, we get the 
asymptotic expansion 

(101) 

The integral then scales as 



(102) 

Here, the spherical integral contributes the factor p, and the 
projection operator contributes p^^. In this limit, the deriva- 
tive 9-y,Xfc does not contribute a factor of p^^, because the 
only factor that can be x^ dependent is the factor that gives 
This factor is dependent on £ and independent of 
p, again because £ is the smallest distance. On the other hand, 
the effect of the derivative d^^^k is to vanish the $2 factor alto- 
gether via an incompressibility cancellation. To see this, note 
that the fusion rule corresponding to the geometry of Fig. |4] 
gives 

Fn+i ~ {wp{xk ~ y,x'fe - y)wy{yik - y,x'fe - y)) $„_i 

(103) 

$2(xfe - y,x'fc - y,Xfc - y,x'fc - y)$„_i, 

(104) 

and from the incompressibility condition we get the tensor 
structure of <1>2 which is 



$2 



(2 + £,,1+1, 2)^ I3j — 61 + 1: 



£jJl 

■ £2 



£?.+ l,2 



(105) 



with £ = ||xfc — x'fc II . The integral 3 can be rewritten as 

J ~ j dy Pafc/3(y)9^,xfc4'2(xfc-y,x'fc-y,Xfc-y,x'fe- 

(106) 

Again, <I>„-i is independent of Xfe and can be pulled out of the 
derivative c^^.x^- However, differentiating with respect to x^ 
wiggles only one of two points (that is Xfe— y, but not x'fc — y), 
which makes it, by chain rule, a derivative with respect to £, 
which in turn vanishes due to the tensor structure of $2 above. 
As a result, we pick the factor cip^^ from the next order term, 
and the integral scales as: 



3 



Cn + l-C>i + l,2-Cn + l,n-l-lg^p-l 



dp p 

Cn+l— Cre + l,2-?n+l,n-l-l 



(107) 
(108) 
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x/c - y 



FIG. 3: UV limit for the case Xfe — s. We employ the fusion rule 
shown in Fig.[T] 



-p-l|y|| 




FIG. 4: IR limit ||y|| oo. We employ the fusion rule shown in 
Fig.S 



the scaling exponent ^„.2 be positive. For a downscale cas- 
cade this gives the condition C2 > and for an upscale cas- 
cade, the condition Cn+i — Cn-i > 0- The assumption of the 
regular fusion scaling is sufficient for that, for both upscale 
and downscale cascades. 

Let us now consider the case where the fusion rules are vi- 
olated according to 

^np = Cp + A^„p (downscale), (113) 

£.np = Cn - Cn-p + Mnp (upSCale). (1 14) 

As we have argued above, as long as the violation is regular, 
UV locally is still maintained. For IR locality, the sufficient 
condition becomes 

A^„+i,2 + A^n+i.n-i > (downscale), (115) 
A^„+i,2 + A^„+i,„_i < (upscale). (116) 

We see that locality survives even the violation of the fusion 
rule hypothesis if A^„+i 2 and A^n+i „_i are both positive 
downscale and negative upscale. 



The locality condition for this integral is Cn+i — — 
^„+i,„_i < and thus Cn+i < Cn+1.2 + Cn+i.n-i- For a 
downscale cascade, the fusion rules hypothesis gives the con- 
dition Cn+i < C2 + Cri-i- For an upscale cascade, the fusion 
rules hypothesis reads = — Cn-p, therefore the condi- 
tion now reads Cn+i > C2 + Ci-i- The condition for locality 
is the same as in the downscale cascade, but the direction of 
the inequality is reversed. 

D. Summary 

Let us now summarize what has been proved. We have 
shown that for either a downscale or an upscale cascade the 
locality conditions are 

(UV): ^„+i,2 > 0, G N, n > 1 (109) 
(IR): Cn+i < en+1,2 + e«+i,n-i, Vn £ N,n > 1. (110) 

for UV locality and IR locality correspondingly. For a down- 
scale cascade, the IR locality condition is satisfied under the 
fusion rules hypothesis 

Cnp = Cp, yp,neN,n>l,2<p<n-2, (111) 

due to the Holder inequality Cn+i < C2 + Cn-i for the scal- 
ing exponents Cn- For an upscale cascade, the fusion rules 
hypothesis gives 

^np = Cn - Cn-p, Vp, n G N, n > 1, 2 < p < n - 2, (112) 

and the IR locality condition is reduced to Cn > C2 + Cn-2 
which is still satisfied, because the Holder inequality reverses 
its direction when the cascade is upscale (see appendix iDb. 
The UV locality condition is also satisfied, but does not re- 
quire the fusion rules hypothesis. All that is required is that 



V. STABILITY OF THE UPSCALE AND DOWNSCALE 
CASCADE 

We now turn to the question of statistical stability with 
respect to forcing perturbations. Statistical stability is de- 
fined as the requirement that there should be a region 3n 
such that (3„({X}„) is negligible relative to contributions to 
Dkn{{^}n) for all {X}„ e 3n in that region. Even when 
the forcing spectrum is confined to a narrow range of scales, 
it is not self-evident that this requirement is satisfied, due to 
feedback loops of F„ onto Qn (see below). 

The first explicit proof that the inertial range of three- 
dimensional turbulence is statistically stable was given by 
L'vov and Procaccia in section II-C-3 of lf52ll . The proof used 
the balance equations of the standard structure functions (not 
the generalized structure functions used in this paper), and it 
covers the case of stability with respect to gaussian forcing 
when the scaling exponents Cn take Kolmogorov scaling val- 
ues Cn = n/i. The value of this proof has gone by unnoticed 
because experiments and numerical simulations have estab- 
lished the statistical stability of the three-dimensional energy 
range beyond all doubt. For the problem of two-dimensional 
turbulence however, where the lack of robustness of the up- 
scale and downscale cascades is the unresolved problem, the 
method used by L'vov and Procaccia in that proof is very il- 
luminating. The main idea is to estimate the scaling exponent 
of the ratio Qn/Dkn and require the appropriate constraint on 
that exponent such that the ratio vanishes asymptotically in 
the inertial range, in the limit of extending the range. 

Consider a geometry of velocity differences {x}„ such that 
all point to point distances have order of magnitude 1, and 
define the scaling exponent (/„ by 

Q„(i?) = g„(i?{x}„) - (^^y . (117) 
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with R a scale in the inertial range. From locality (proved in 
the previous section) we also know that. 



Dkn{R) = Dkn{R{^}n) 



It follows that the ratio Qn/Dkn scales as 



Dkn{R) 



R 



<?„-(a+i-i) 



(118) 



(119) 



In a direct cascade, such as the energy cascade of three- 
dimensional turbulence and the enstrophy cascade of two- 
dimensional turbulence, this ratio must vanish in the limit 
^0 +00. It follows that the condition for the statistical 
stability of a downscale cascade reads 



(C„+i - 1) > 0, Vn G N, n > 1. (120) 



In an upscale cascade, such as the inverse energy cascade of 
two-dimensional turbulence, the same ratio must vanish in the 
limit ^0^0. This leads to the same condition with the in- 
equality reversed: 

Ag„ = (z„ - (C„+i - 1) < 0, Vn e N, n > 1. (121) 



A. The case of gaussian forcing 

For the simplest case of Gaussian delta-correlated in time 
forcing, the exponents Qn can be calculated in terms of (n- 
This makes it possible to investigate statistical stability rigor- 
ously. 

We begin with the assumption that is a delta-correlated 
stationary gaussian field with (/^(x)) = 0, and 

(/a(xi,ii)/^(x2,i2)) = 2eCa/3(xi,X2)5(tl -t2), (122) 



where e is constant, and Cq.^ is normalized such that 
Cq.q(x, x) ~ 1. Without loss of generahty we may assume 
that da fa = 0, and therefore 'S'apfp = Thus, we have the 
identity 

j (iy ^'/37(x2 - y)CQ^(xi,y) Cq/3(xi,X2), (123) 

which will be used below. 

We define the forcing scale £o from the Taylor expansion 



C.Mx + y,x).^-4^)(M') +0(0, (124) 



valid in the limit ||y|| <C £o- Note that the odd-order terms 
vanish by incremental isotropy. In the limit ||y|| ^ i^, on the 
other hand, we have the asymptotic expansion 



Ca/3(x+y,x) - ( — 



(125) 

Note that a, which is an unspecified scaling exponent depen- 
dent on our choice of stochastic forcing, must satisfy a > 0, 
since the correlation must vanish at ||y|| +00. Also 
note that e is the total rate of energy injection. In general, 
the work done on the fluid is £i,i(x) = fa{x)ua{x). For 
delta-correlated forcing, it is easy to show that (ej„(x)) = 
eCaQ,(x, x) (see proof in appendixICb. 

Recall that the total forcing term Q„ is given by 



Q„({X}„, t)=J2 Qknimt Xfc, t), (126) 
where Qkn reads 



fc=i 



Qr^-""-^'({x}„-i,Y,t) 



n w„,(Xfe,t) 



(127) 



with (/5q(X, t) = /q(x, t) — fa{x', t). For Gaussian forcing, it can be shown (see appendix iBt that the forcing contributions 
Qkn to the generahzed balance equations read 



Tt-l 



(128) 



1=1 

with Qq/3(X, Y) given by 

Q„;3(X,Y) = {wa(X,t)ipp{Y,t)) ^2e y dz [P„^(x - z) - P„^(x' - z)][C/3^(y, z) - C/3^(y', z)] (129) 



2e[Capiy,^) - Capiv',^) - a/3(y,x') + Cat3{y',^% 

I 



(130) 



The physical intuition is that there is a feedback loop between forcing, whose spectrum is defined by Qap (X, Y), and the re- 
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suiting behavior of turbulence which is captured by the struc- 
ture functions F„. More specifically, we see that Fn~2 pro- 
vides feedback to Qn, when the forcing is gaussian. For sta- 
tistical stability we need this feedback to be negligible in the 
inertial range. 

The immediate implication of eq. (11281 ) is that = Cn-2 + 
52 with 52 the scaling exponent of Qap- It follows that 

Ag„ = (C„-2 + 92) - (Cn+i - 1). (131) 

The remaining challenge is to calculate (72. We will see that 
(72 depends on whether the cascade is upscale or downscale. 



which suggests that for a downscale cascade, q2 = 2. It is 
easy to see that for a monofractal velocity field with (n = nh, 
the stabihty condition reads 

Ag„ = (C„-2 + 2) - (C„+i - 1) (136) 
= 3- 3/i > 0,V7i e N : n > 1, (137) 

which requires h < 1. In a multifractal case one has a linear 
combination of independent monofractal contributions, and it 
can be shown that the constraint < < 3 is a sufficient 
condition for statistical stability. This follows from the in- 
equality Cri+i < Ca + Cri-2 (see appendix iDb: 

Ag„ = C„-2 - Cn+i + 3 (138) 

> Cn-2 - Cn-2 - Ca + 3 (139) 

= 3-C3 > 0,V7i e N : 71 > 2. (140) 

For n = 2, we get Aq2 = 92 ~ (Cs ~ 1) = 3 — (3, which 
implies, from the stability condition Aq2 > 0, that < Ca < 3 
is also a necessary condition. 

For the case of the downscale energy cascade of three- 
dimensional turbulence we have Cs = 1 , which can be derived 
from the solvability condition for the homogeneous equation 
02^3 = imSllil]. This satisfies the sufficient condi- 
tion < Cs < 3 for statistical stability very generously, so 
it is hardly a surprise that the energy cascade is so robust. 
Also worth noting is that for a hypothetical downscale helic- 
ity cascade we have Cs ~ 2 , which also satisfies the stability 
condition. 

The story changes for the case of the downscale enstrophy 
cascade. We know, from combining the Eyink and FaUcovich- 
Lebedev theories of the two-dimensional enstrophy cascade 
ll22l I23I F74I F75I1 . that when it exists with constant enstrophy 
flux, the enstrophy cascade has no intermittency corrections. 
Thus, the scaling exponents Cn all satisfy the monofractal 



In the rest of this section, we will derive the separate stability 
conditions for a downscale cascade and for an upscale cas- 
cade. 



B. Stability conditions for downscale cascades 

For the case of a downscale cascade, using the Taylor ex- 
pansion of (5q,^(X, Y) in the limit ||X — Y|| 0, the scaling 
of Qap can be estimated as 



(132) 
(133) 
(134) 
(135) 



scaling Cn = , which implies that 

Ag„ = Cn-2 - Cn+1 +3 = 0. (141) 

So, we have a borderline situation where the stability condi- 
tion is neither satisfied nor broken! Consequently, the actual 
stability of the downscale enstrophy cascade is not decided by 
scaling exponents but by the numerical coefficients in front of 
Qn and Dkn ■ This is where it gets interesting. 

The leading contribution to Qn is proportional to the to- 
tal rate of energy injection e. However, one should bear in 
mind that the downscale enstrophy cascade is forced by the 
combined effect of both the forcing term and the large- 
scale dissipation term (— 1)'"+^/3V^^'"mq. As a result of this 
combined forcing, the enstrophy cascade is injected with a 
smaller enstrophy rate rjuv and a very small energy rate e„„ 
with rjuv < V and Suv ^ £■ If we assume that this com- 
bined effect itself can be modelled as gaussian forcing, then 
the leading contribution to the effective forcing on the enstro- 
phy cascade is proportional only to the rate Suv of the sublead- 
ing downscale energy flux. Because e„„ vanishes rapidly as 
the separation of scales in the enstrophy cascade is increased 
ll36ll40ll . this leading contribution can be made as small as de- 
sired simply by taking the limit ^ 0+. For small enough 
downscale energy flux the next order term with 92 > 3 
becomes dominant, and combined with Cn = it is easy to 
show that the stability condition is now Aqn > 0. 

The conclusion from this analysis is that the stability of the 
downscale enstrophy cascade requires that the accompanying 
downscale energy flux should be very small. For that to hap- 
pen, we need two things: First, it is necessary to have a dis- 
sipation sink at large scales to absorb most of the injected en- 
ergy at the forcing scale or at larger scales. Second, we must 
have a large separation of scales between the forcing scale 
and the dissipation scale at small scales, which means that a 
significant amount of numerical resolution is required. These 



Q„/5(X, Y) = 2£[C„^(y,x) - C„/3(y',x) - Ca^Y, x') + C^My', x')] 

= 2£[(C„^(y,x) - C„^(x,x)) - (C„;3(y',x) - C^pi^,^)) - (C^piy,^') - C„^(x',x')) 

+ (C„;3(y',x')-aMx',x'))] 
~ {2s/il)[\\y - xf - ||y' - xf - ||y - x'f + ||y' - x'f] ^ 
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two requirements, we believe, are the reason why it has been 
so difficult to reproduce the enstrophy cascade in numerical 
simulations. It is worth noting that Tran and Bowman ll2lll 
came to a similar conclusion by a different argument, that the 
robustness of the downscale enstrophy cascade requires a van- 
ishing downscale energy flux. 

C. Stability conditions for upscale cascades 

The fundamental difference between an upscale cascade 
and a downscale cascade with respect to stability is that in 



which gives (72 = —a < 0. For a monofractal velocity field 
with (n = nh, the stability condition reads 

Aqn = Qn - (Cn+1 - 1) (144) 
= 92 + Cn-2 - (Cn+1 - 1) (145) 

= (72 + l-3/i<0,VneN:n> 1. (146) 

Since q2 < 0, the condition h > 1/3 is sufficient. For the 
more general multifractal case, using the inequality Cn+i > 
Cs + Cn-2i we can upper-bound Aq^ as follows; 

Aqn =q2+ Cn-2 ^ (Cn+1 - 1) (147) 
< 92 + 1- (Cn-2 +C3)+ Cn-2 (148) 

= 92 + 1 -C3,Vn e N : n > 2. (149) 

For n = 2, we get an equality: Aq2 = 92 — (Cs — 1)- Thus, 
the stability condition Aqn < is satisfied when Cs > 92 + 
1, which is indeed satisfied when (j2 < and Cs > 1- Fot 
the inverse energy cascade of two-dimensional turbulence, we 
have C3 = 1 which satisfies the requirements for stability. 

There is however another effect that can destabilize the in- 
verse energy cascade. We have shown in a previous paper ll47ll 
that the loss of asymptotic homogeneity by the effect of the 
boundary conditions on the flow amplifies the sweeping term 
/„ at the large scales. As a result, at sufficiently large length 
scales, the ratio In/Dkn becomes significant, and excites a 
particular solution superimposed on top of the homogeneous 
solution associated with the inverse cascade. The particu- 
lar solution corresponds to the coherent structures associated 
with the "energy condensation effect". The formation of these 
coherent structures is very likely to further intensify the ratio 
In/Dkn- As we have explained in the introduction, it has been 
shown that if these coherent vortices are removed before the 
evaluation of the energy spectrum, the usual inverse energy 
cascade spectrum is recovered ll30l[32i[34ll . This result is con- 
sistent with our theory, and it confirms that the homogeneous 



the upscale cascade the Holder inequalities now take the form 
Cn+fc ^ Cn + Cfc' ^nd the condition for Statistical Stability reads 
Aqn < 0, Vn > 1. We will now prove that inverse cascades 
are always statistically stable with respect to variations in the 
forcing statistics, provided that Ca > 1- This is consistent 
with the numerical evidence that the inverse energy cascade is 
much easier to obtain in simulations than the direct enstrophy 
cascade. 



Again, using Taylor expansion in the limit ||X — Y|| 00, 
we see that Qap scales as 

(142) 
(143) 

I 

solution, corresponding to the inverse energy cascade, exists 
side by side with the particular solution, corresponding to the 
coherent structures, even when the particular solution is dom- 
inant. The possible role of the sweeping term on the stability 
of the enstrophy cascade is currently not well-understood. 

VI. CONCLUSION AND DISCUSSION 

We have shown that the non-perturbative locality of the in- 
ertial ranges of two-dimensional turbulence is an immediate 
consequence of the fusion rules hypothesis. The physical in- 
terpretation of what we have done is to prove, strictly in the 
context of the incompressible Navier-Stokes equation, that 
universality implies locality. A proof of the fusion rules by 
diagrammatic theory is essentially the converse and more in- 
teresting claim: that locality implies universality. This result 
leads to an apparently curious paradox: the usual understand- 
ing of locality, in terms of triad interactions in Fourier space, 
suggests that a necessary condition for locality is that the en- 
ergy spectrum E(k) must have slope between and k~^. 
This corresponds to the inequality < C2 < 2. The paradox 
is that this constraint does not appear anywhere in our locality 
proof! In recent work, Eyink IITql investigated the locality of 
the downscale enstrophy cascade and the inverse energy cas- 
cade using a filtering method IttI FtsI [79ll . His argument also 
leads to the inequality < C2 < 2 as a sufficient locality con- 
dition. It follows that whereas the inverse energy cascade is 
local the direct enstrophy cascade is IR marginally-nonlocal. 
Unlike the argument in this paper, Eyink's argument has only 
considered the kinematic locality of the flux term and not the 
statistical locality associated with unfused higher-order struc- 
ture functions. On the other hand, our argument is less rigor- 
ous in its present form, as it assumes the fusion rules without 
proof. 

A fundamental problem with establishing locality in 



)„;3(X,Y) = 2e[C„;3(y,x) -Ca;3(y',x) -C„;3(y,x') + a/3(y',x')] 



2e 



|y - x| 



|y' -x| 



|y - X' 



|y' - X' 
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Fourier space is that the Fourier transform involves an inte- 
gral that ranges over every length scale, including the forcing 
length scales and the dissipation length scales. To preserve 
locality, the main contribution to the integral must come from 
the inertial range. The inequality < (2 < 2 comes in as 
a necessary condition for the survival of locality under the 
Fourier integral |63l- The same issue arises when locality is 
characterized with a filtering transform (i.e. forward Fourier, 
truncation, backward Fourier), as was done by Eyink llTdl , al- 
beit with a broader definition of filtering. Beyond that, the un- 
derlying argument based on diagrammatic theory ifsoHsTl I52I1 
that justifies the fusion rules hypothesis itself can impose fur- 
ther constraints on ^2, which still need to be investigated care- 
fully. For example, one other way the constraint < C2 < 2 
can come in is if we require perturbative locality for each 
Feynman diagram ifsoll . Perturbative locality may be a nec- 
essary condition for the fusion rules hypothesis. If that is true, 
then perhaps < (^2 < 2 is implicitly assumed when we pos- 
tulate the fusion rules hypothesis. We have also shown in this 
paper that the related condition < Cs < 3 is required for sta- 
bility under Gaussian forcing, which is as essential as locality 
for the existence of a universal inertial range. 

It should be stressed that any constraints on scaling ex- 
ponents needed only to prove the fusion rules hypothesis by 
Feynman diagrams , are needed only to establish the univer- 
sality of the scaling exponents C„ of the inertial range. We 
should expect to find that the conditions for locality itself are 
weaker For example, we have shown in this paper that local- 
ity is possible even when the fusion rules fail, provided that 
the fusion exponents ^„p deviate in the correct direction. In 
fact, it is possible to have local interactions, as per our def- 
inition, even when the underlying diagrammatic theory does 
not yield local Feynman diagrams! This scenario is not en- 
tirely hypothetical; in the case studied by Refill of an en- 
strophy range under strong Ekman dissipation, this may be 
precisely what happens, with the slope being steeper than k^'^ 
and non-universal, but still allowing an appearently local en- 
strophy cascade to exist. 

The key idea that can help us unravel these paradoxes is that 
the non-perturbative locaUty studied in this paper is a weaker 
condition than perturbative locality. Non-perturbative locality 
requires only the combined effect of all Feynman diagrams to 
be local. Perturbative locality, on the other hand, requires that 
each diagram individually should be local. This distinction 
between perturbative and non-perturbative locality may clar- 
ify the paradoxical situation with the enstrophy cascade where 
the spectrum of the enstrophy cascade is consistent with a di- 
mensional analysis argument based on a locality assumption 
even though the slope is too steep to be self-consistent with 
that assumption! Adding a logarithmic correction resolves the 
situation in a one-loop closure model Isill . and the combina- 
tion of more recent results by Falkovich and Lebedev 1I23I1 and 
Eyink [75] suggest that the same logarithmic correction per- 
sists for the exact theory, with no higher-order adjustments. 
Nevertheless, a reconcilliation of the spectrum slope and the 
locality requirement is still an "uncomfortable" notion, to say 
the least. We believe that a possible resolution of this para- 
dox is to claim that the enstrophy cascade is local in the non- 



perturbative sense, as far as the exact theory is concerned, and 
borderline non-local only in the perturbative sense. From a 
physical standpoint the relevant locality needed as a precon- 
dition for establishing the existence of an inertial range is the 
non-perturbative locality. However, some confusion can arise 
from the fact that closure models unwittingly exchange non- 
perturbative locality with perturbative locaUty! 

The careful reader will note that the non-perturbative lo- 
cality is also weaker than the more intuitive (and less rigor- 
ous) physical understanding of locality as the notion that the 
effect of the forcing range and dissipation range is "forgot- 
ten" in the inertial range. We may designate locality, in this 
sense, as "strong" locality, so that it can be distinguished from 
the weaker non-perturbative statistical locality. The proposed 
theory can help make the meaning of this notion of "strong" 
locality more rigorous. The key idea is that it is possible 
to have local interaction integrals in the contributions to the 
0„F„+i term of the balance equations and still pick up an ef- 
fect from the forcing range or the dissipation range into the 
multi-dimensional regions dn that are supposed to be the in- 
ertial range, in our generalized sense. It all depends on how 
much forcing and dissipation "wish to creep into" the iner- 
tial range. We can find that out by comparing the magnitude 
of the Qn , In, and DnFn terms of the generalized balance 
equations against the magnitude of the contributions Dkn to 
the interaction term. Thus, we find that there are three distinct 
conditions that need to hold to have strong locality: first, the 
interaction integral itself has to be local; second, we need to 
establish the property of statistical stability which will guar- 
antee that the forcing effect Qn and the sweeping interactions 
/„ do not creep into the inertial range; third, a calculation of 
the shape of the dissipation range can show whether there is 
a wide enough region 3n in which the dissipation term DnFn 
is negligible. One advantage of the generalized balance equa- 
tions framework is that it allows us to account mathematically 
for these three distinct effects separately. 

In this paper, we examined only the first condition and part 
of the second condition. We have shown that statistical sta- 
bility with respect to forcing applies unconditionally for the 
inverse energy cascade. For the enstrophy cascade, statisti- 
cal stability requires large-scale dissipation and a vanishing 
downscale energy dissipation. For any downscale cascade in 
general, stability constrains the corresponding Holder expo- 
nent as h < 1. For an upscale cascade, the corresponding con- 
straint is h > 1 /3. We began considering the role of sweeping 
in a previous paper ll47ll . and the role of the dissipation term 
will be studied in future work. 
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APPENDIX A: DERIVATION OF THE BALANCE 
EQUATIONS 

In this appendix we give a detailed derivation of the gener- 
alized balance equations. Recall that we defined the general- 
ized structure function F„ as 



i^„({X}„,t) = 



Lfe=l 



(Al) 



By differentiating Fn with respect to t and substituting the Navier-Stokes equations we obtain: 

n 



dFn{t) _ ^ / 5wQfc(Xfe,x'fe,t) 



dt 



fc=i 



dt 



(A2) 



fe=i 



Here, the terms z/ J„ and are the contributions of the small-scale and large-scale sinks with 

n 

j:^"--""({x,x'}„,t) =^(V^^: + v2«jF„({x,x'}„,i), 

n 

ij:^"--""({x,x'}„,o = ^(v^f-- + v^:'pF,,i{^y}„,t), 



(A3) 
(A4) 



fc=i 



where V^^ is the Laplacian with respect to x^; V^/^ is the Laplacian with respect to x'^. Also, Nkn represents the contributions 
of J'q,^9^(u^w^), and Qkn represents the contributions of Vapfp, and they read: 



o^r'""({x,x'}„,t) 



7V,"r^-""({x,x'}„,t) 



n 

Jl Wa,(xi,x'i,i) 

n 

Y[ WQ,(xi,x'i,i) 

n 

Y[ Wa,(Xi,x'i,i) 

l=l,l^k 



yaM{^k,t)~M^'k,t))j, 

5'afc/3[9T.,xt,(M/3,XfcM7,xJ " 9^,x', (w/3,x'fc "7,x'fc 
5'afc/3^/3(Xfc,x'fc,t)^ 

y dyPa,/3(y)^/3(xfc -y,x'A,. -y,t) V 



(A5) 



(A6) 



(A7) 



(A8) 



Here we use the abbreviations Ua,Xf, = Ua{xk,t) and Wq.x'j. = Ua{x'k,t), Wa,k = '"^^(xfc, x'^, t), and da.Xk is the spatial 
derivative in the a direction with respect to x^. Also, ^^^(xfc, x'^, t) is the non-linear factor defined as: 

J^p{xk,X.'k,t) = (97,xfc(w,3,xfc%,xj - 9^,x'fc(lt/3,x'fcU7,x'J = U7,XA,<9^,Xfc(u/3,Xfc - lt/3,x' J + ^^x'^ 9^,x'fc (M,3,Xfc " W^,x'J 
= S^.Xfcl^.XfcW^^fe) +9^,x'fc(tt7,x'fcW/3,fe) = %,x,9^,XfcW^,fe +%,x'fc<9^,x'fcW/3,fc■ 



It is easy to see that the nonlinear terms Nkn cannot be writ- 
ten exclusively in terms of velocity differences. The remark- 
able characteristic of the derivation of the balance equations 
by L'vov and Procaccia ll43ll is that the nonlinear term Nkn 
is rearranged as the sum of a local term Dkn and a sweep- 
ing term I^n such that the local term can be expressed as a 



r 



linear operator on Fn+i- Although L'vov and Procaccia II43I1 
eliminated the sweeping term on the grounds of global ho- 
mogeneity, we believe it is appropriate to retain it here in its 
simplified form. 
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To isolate the sweeping term, we define a generalized mean 
velocity Uq({z, z'}„, t) as: 

1 " 

Ua{{z,z'}n,t) = — VK(Zfc,i) +M„(z'fc,i)), (A9) 
/n. ^ — ^ 



and the corresponding velocity fluctuation 

1 " 

Wa(x, {z,z'}„,t) =UQ(x,t) - 'Ua({z,z'}„,t) = — ^ [iOq (x, Zfc) + (x, z'j 



2?! 

We may then decompose Kct(xfc, x'^, t), in general, to 

]NfQ(xfc,x'fc,t) = §a(xfc,x'fc, {z,z'}„,t) +ila(xfc,x'fc,{z, z'}„,t), 

where Sq. and are defined as: 

§Q(xfc,x'fc,{z, Z'}„,t) = lX/3({z,z'}„,t)(9^.x^ + 9/3,x'JwQ(xfc,x'fc,t), 

£Q(xfc,x'fc,{z, z'}„,t) = [w/3(xfc, {z,z'}„,t)5/3,xfc + w/3(x'fc, {z,z'}„,t)93^x'JwQ(xfc,x'fc,i). 
In general, {z, z'}„ can be chosen any way we wish. Here, we specifically use the choice: 

?\f^(xfe,x'fe,t) = iL^(xfc,x'fc, {x,x'}„,t) + S^(xfe,x'fe, {x,x'}„,i). 
This gives the decomposition N^^n = Dkn + Ikn with 



(AlO) 



(All) 



(A12) 



(A13) 



i^^r-""({x,x}„,t) 



/^r^-""(Kx}„,t) 



Wq, (x;,x'(,i) 

n 

J]^ Wq, (Xi,x'i,t) 

l = l,l^k 



j 'iy^'Qfc/3(y)''^/3(xfc - y,x'fc - y, {x,x'}„,t)^ , (A14) 
y" (iyPafc/3(y)§/3(xfe -y,x'fc -y, {x,x'}„,t)\ . 



(A15) 



Here /fe„ represents the sweeping interactions and Dkn represents the local interactions. 



The sweeping term Ikn can be simplified as follows: We use the decomposition Pq,^(x) = (5q^(5(x) — P^^(x) to split I^n to 

two terms: lun = lin + ^hn ^^'-^ ^hn corresponding to SapS{x) and ij^J corresponding to P^^(x). We also use J'lj.^wa = 0. 
The integral inside the ensemble average of Ikn splits to two parts: Ii and l2- The first part Ii reads: 

h= j dy 6a^p5{y)§p{xk - y,x'fc - y, {x,x'}„,t) = §ak{^k,^'k, {x,x'}„,t) 



= lX^({x,x'}„,t)(5^,xfc + 5^,x'Jit;Qfc(xfc,x'fc,t). 
The second part I2 is shown to be zero by incompressibility: 



h = J dyPlf^{y)§p{xk -y.^'k - y, {x, x'}„, 

y (iyP'j^(y)l(^({x,x'}„,t)(5^,Xfc + 9^,x' Ju;/3(xfc -y,x'fc -y,t) 
= 'U^({x,x'}„,t)(a^,xfc +3^,x'J / dyPf^^(y)w/3(xfc -y,x'A; -y,t) =0. 



(A16) 
(A17) 

(A18) 
(A19) 
(A20) 



Because is the nonlocal part of the projection operator, this result implies that the pressure effect does not contribute to the 
sweeping interactions or to the violation of incremental homogeneity. Thus, Ikn is determined by /i and it simplifies to 



/^r^-""({x,x'}„,t) = 



l=l,l^k 



l[^({x,x'}„,i)(a^_xfc + S^.xvJifQjxfejx'fejt) 



(5^,xfc +a^,x'J ( ii^({x,x'}„,t) 



.1=1 



(A21) 
(A22) 
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This result was given previously by L'vov and Procaccia in section IV-B and appendix B of Refl43[ 

We will now show that the local interaction term Dkn can be written as a linear transformation of Fn+i- First, note that 

£c<(xfc,x'fc,{x,x'}„,i) = [w^(xfe, {x,x'}„,t)5/3,xfc +W/3(x'fe,{x,x'}„,t)5;3,x'fc]wa(xfc,x'fc,i) (A23) 

= 9/3,x& [w/3(Xfe, {X, x'}„, t)Wa{:>Ck,^'k,t)] + dp^^'^ [w/3(x'fc, {X, x'}„, t)Wa{Xk, x'fc, t)] (A24) 



1 " 

= — ^9/3,xJ(w/3(xfc,Xi,t) +u7^(xfe,x';,i))u;Q(xfc,x'fe,t)] (A25) 
1=1 

1 " 

+ ^^9/3,x'fc[(w/3(x'fc,Xi,t) + u;^(x'fc,x';,i))wQ(xfc,x'fc,t)], (A26) 
1=1 

which gives: 

1 " 

La{xk -y,x'fc -y,{x,x'}„,t) = -^^dfi^^^iiwpixk -y,x;,i) +u;/3(xfc - y, x',, t))u;a(xfc -y,x'fc -y,t)] (A27) 

1=1 

1 " 

+ ■^^9p^^>^[{wi3{-x.'k -y,x;,<) +u;,3(x'fc - y, x';, t))wa(xfc -y,x'A; -y,t)]. 
1=1 

(A28) 

It follows from substituting the above to (IA14I ) that Dkn is given by 

i?^r -""({x,x'}„,t) - f E / dyP.,My)^."r ' "'"-^'' ' ""({x,x'}„,y,t), (A29) 

" 1=1 •' 

with Dkni = Dknii + Dkni2 + Dknia + Dknii, and 

^fcnZi (ix,x|„,y, ^) - o'a„_n,xfc^„+i (|X„i)„^;^,Xfc - y,x - y, iX„j),„^fc+i,Xfc - y,x;j, 

(A30) 

^fc„i2 (ix,x |„,y,tj - o'a„_^i,xfci'„+i (|X™),„^;^,Xfc - y,x fe - y, |X™)„^j.+i, x^ - y,x 

(A31) 

-'^fcnia llX, X |„,y,tj — Ca„_^i,x'fc^„+i a^mj-,„=i,Xi. - y,X - y, |J^TOj„j^j,^;^,X k - y,x;j, 

(A32) 

J->knU (ix,x |„,y,tj - Ca„+i,x'fe-f'„+i (|X,„)„^;^,Xfc - y,x fc - y, |X,„)„^j,+i, x ^ - y,x ;) 

(A33) 



APPENDIX B: FORCING CONTRIBUTION FOR GAUSSIAN exploit the following mathematical result: if /„ (xi , ti ) is a 
FORCING Gaussian stochastic field, the ensemble averages of the form 

(/a(xi, ti)R[f]) can be evaluated for any analytic functional 
We give here a proof of equations ( fT28] ) and ( flJOl l. closely R[f ] by the following integral 
following the argument in section II-C-3 of Ref. El We 



(/„(xi,ti)i?[/]) = j d^2dt2 (/o(xi,ii)/^(x2,i2)) 

We begin the proof by defining the following response functions 



SR[f] 



ai ■■■am/3i---/3„ 



nai- 



Sap(X.,ti;y,t2) 



({X}„,t, {y,r}„) 



SMy,t2) 



n 



(Bl) 

(B2) 
(B3) 



.1=1 
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For the case ii = t2 = t, the response function Sq^(X, t; y, t) is given by 

g„;3(X, t; y, t) = (l/2)[F„;3(x - y) - P„/3(x' - y)]. 



(B4) 



This is proved in appendixlcl Likewise, for the case m — \ and ri = t, the response function S"^ ""'^({X}„, t, y, t) is given 
by 



S?--""''({X}„,i,y,t) = 



\^//3(y,i) 



i=i 



E 

fc=i 



(5Wa^^(Xfc,t) 



(B5) 
(B6) 



k=l 



Here we exploit the fact, first pointed out in Ref. 52, that the variational derivative {Swa^ (X^., <))/((5/^(y, t)) is not correlated 
with the velocity differences Wq, (X; , t) because no time is being allowed for interaction to develop a correlation. Using (IB II) 
the correlation between Wa (X) and //3(y) is given by 



{fl3{y,t)fyiz,T)) ^2e / dz9a'r(X,t;z,t)Cfiy{y,z), 



and it follows that 



g„^(X,Y)=2£ J dz9a^{X,t;z,t)[Cp^{y,z) - Cp^{y',z)] 

= 2eJ dz[Pajix-z)- Pc^{x' -z)][Cp^iy,z)-Cpjiy',z)] 
Using a similar argument for the more general case, we get 

l[wa,iXi,t) My,t)\ ^ Jdz Jdr (jj^^ n^-'(^''^) 



_<5/^(z,r) 

2e J dz g^-:'^"-^'^({X}„_i,z)C3^(y,z), 



(/a(y,t)/^(z,r)) 



(B7) 

(B8) 
(B9) 

(BIO) 

(BID 
(B12) 



and it follows that 



1=1 



(B13) 
(B14) 



{My,t)~My',t))j 

2e J dzg^-:'^"-''({X}„_i,z)[C^,(y,z)-C^,(y',z)] 

n — l „ 

2eE^"-i2'"'"'"'^""""-'({X}Li) / dz<5^,,[Xut-zMCp,[y,z)~Cp,{y',z)] (B15) 

1=1 

n-l 

E^"2'"'-^"'""''""-^({X}Ll)Qa,MX/,Y). 



(B16) 



1=1 



This concludes the proof. 



APPENDIX C: EVALUATION OF THE ONE-TIME 
RESPONSE FUNCTION 



We show how to calculate the one-time response function 
and use it to show that the ensemble average of the rate of 



r 



energy injection £„i(x) is given by (ei„ (x)) = £Caa(x,x). 
This argument was given previously by McComb ll82ll . 
We begin with the definition of the response function: 



LTa/3(Xi,ti;x2,r2j — ( YTl TT 

\ 0//3(X2,l2j 



(CI) 
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We first show that at equal times ti = t2. Gap is given by 



To show this, note that from linearity with respect to forcing 



GQ/3(xi,<i;x2,ti) = iPa/3(xi -X2). (C2) 



r* /■* 

MQ(x,t) = Ua(x, 0) + / ds Aa[ua{s)]{r) + ds I dy Paf3{x.~ y)fa{y,s), 











(C3) 



where Aa[ua{s)]{r) represents the effect of the advection and pressure term. For convenience, we use the abbreviation 
^q/j/q- It follows that 





[ ' ds 


SAa [■tiQ(s)](xi) 


Sfi3(x2,t2) 


Jo 






[ ' ds 


SAa[ua{s)]{r) ^ 




Jt2 


Sfp{x2,t2) 




[ ' ds 


6Aa[ua{s)]{r) 




Jt2 


Sfl3{^2,t2) 



ds dy Pa(3{yi~y)fp{y,s) 



(C4) 



Xf / . ^ / / dt[Hih-t}Papi^i-y)]f,3iy,t) (C5) 



Hih^t2)Papixi~X2)- 
I 



(C6) 



with H{t) the Heaviside function, defined as the integral of a 
delta function; 



H{t) 



t I 1, if t e (0, +00) 

S{t) dr = \ 1/2, if t = . (C7) 

0, ifte(-oo,0) 



For ti ~ t2, the integral of the first term vanishes and H{Q) = 
1/2, therefore it follows that 

GQ/3(xi,ii;x2,<i) = ^Pap{xi -X2). (C8) 

Also note that in fact there is a discontinuity in the response 
function and 

lim Ga/3(xi,t + At;x2,t) = P„^(xi -X2). (C9) 
From this result, it immediately follows that: 

(ez„(x))= J rfxo J dto {fa{x,t)fp{xo,to)) (CIO) 

= J dxo 2eCa/3(x,xo)Ga/3(x,t;xo,t) (Cll) 

= J dxo eCai3{x,xo)Pap{x - xo) = eCaa{x,x). 

(C12) 

APPENDIX D: SCALING EXPONENT INEQUALITIES 

We will show here that for an downscale and upscale cas- 
cade, correspondingly, the scaling exponents satisfy the in- 
equalities 



Cn+k < Cn + Ck (downscale), 
C"+fc > Cn + Cfe (upscale). 



(Dl) 
(D2) 



The first of these inequalities is well-known. The key result 
here is the second inequality, corresponding to the case of an 
upscale cascade, whose direction reverses, thus giving a con- 
vex upward (or flat) dependence of (n as a function of n. This 
should be contrasted with the case of a downscale cascade 
where the dependence of C„ on n is convex downward (or flat). 
The proof is "folklore" and it uses the Schwarz and Holder 
inequalities. An earlier version of this argument was given by 
Frisch ll63il64ll . who in turn cites Feller Issll . 

Letp, q € {1, +00) with 1/p+l/q = 1, and let cj), be two 
random variables with (f> > and ip > 0. The Holder inequal- 
ity for ensemble averages states that {4>ip) < {<jf)^^^ (-0'^)^^^. 
For p = q = 1/2 it reduces to the Schwarz inequality: 

We begin by defining w{R) as the absolute value of the 
longitudinal velocity difference: 



w{R) = \ (u(x + Re, t) - u(x, t)) ■ e\ 



(D3) 



where x G M"^ is given and e is a unit vector The proof is 
based on the following two assumptions: (a) For a downscale 
cascade, in the limit £0 ^ oo- w{R) scales as ([w(i?)]") ~ 



For an upscale cascade, the same scaling law holds 



for the limit £q O"*". (b) For finite io there is a range of 
scales where the above scaling law continues to hold as an 
intermediate asymptotic 

The proof uses two "helper" inequalities that are interest- 
ing in themselves. The first "helper" inequality is deduced by 
choosing (j) = ^, ^ [i(;(i?)]("+i)/2 and 

employing the Schwarz inequality. It follows that 



{HR)r) 



{[wiRT-') {[wiR)r+') . 



(D4) 
(D5) 
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and therefore 



By similar reasoning, we find that 



{[w{RT 



2C„-Cr.-l-C^ + l 



< 1. 



(D6) 

To satisfy this inequality under the limit £o ^ (x we require 
2C?i — Cn-i ^ Cn+1 > 0. Thus we get for a downscale cascade: 



Cri+l - Cri < Cri - Cn-1 (doWnSCalc). 



(D7) 



Likewise, for an upscale cascade, the inequality must be satis- 
fied in the limit £o 0+, whichrequires 2C„ — Cn-i — Cn+i 1^ 
0. Thus, for an upscale cascade we have 



Cn+l - Cn > Cn - Cri-1 (upSCalc). 



(D8) 



The second "helper" inequality is deduced by choosing ~ 
and ip = ~ 1 and employing the Holder 

inequality with p ~ (n + and g = n + 1. It follows that 

= ([«;(i?)]"+i)"/^"+'\ (DIO) 
which implies that 



{[wiRT) 



C„-(n/(n+l))Cn + l 



< 1. 

(Dll) 



n + 1 

Cn+i < C« (downscale), (D12) 

n 



n + 1 

Cn+i > Cn (upscale). 

n 



(D13) 



Now let us consider the case of a downscale cascade. We 
assume with no loss of generality that n > k (otherwise for 
the following step, one may exchange n and k). Combining 
the inequalities (ID7l i and ( ID 121 ) gives 



Cn+fc-Cn= J2 iCa+1 - Ca) < HCk+1 ~ Ck) (D14) 

(D15) 



<fc( ^a-ai - a, 



Thus we estabUsh that 



Cn+fc < Cn + Ck (downscale). 



(D16) 



For the case of the upscale cascade, the exact same argument, 
with every inequality reversed, gives 



Cn+fc > Cn + Cfc (upscale). 



(D17) 
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